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ABSTRACT 

Using the 100-m radio telescope at Effelsberg, we mapped a large area around the 
Andromeda Galaxy in the 21-cm line emission of neutral hydrogen to search for high- 
velocity clouds (HVCs) out to large projected distances in excess of 100 kpc. Our 3 a 
Hi mass sensitivity for the warm neutral medium is 8 x 10'* M©. We can confirm the 
existence of a population of HVCs with typical Hi masses of a few times 10^ M0 near 
the disc of M31. However, we did not detect any HVCs beyond a projected distance 
of about 50 kpc from M31, suggesting that HVCs are generally found in proximity of 
large spiral galaxies at typical distances of a few 10 kpc. 

Comparison with CDM-based models and simulations suggests that only a few of 
the detected HVCs could be associated with primordial dark-matter satellites, whereas 
others are most likely the result of tidal stripping. The lack of clouds beyond a pro- 
jected distance of 50 kpc from M31 is also in conflict with the predictions of recent 
CDM structure formation simulations. A possible solution to this problem could be 
ionisation of the HVCs as a result of decreasing pressure of the ambient coronal gas 
at larger distances from M31. A consequence of this scenario would be the presence of 
hundreds of mainly ionised or pure dark- matter satellites around large spiral galaxies 
Hke the Milky Way and M31. 

Key words: intergalactic medium ~ galaxies: evolution - galaxies: individual: M31. 



1 INTRODUCTION 

One of the currently most favoured cosmological models is 
the so-called Lambda Cold Dark Matter (ACDM) model. 
It assumes that the evolution of the universe is dominated 
by dark energy and dark matter which in sum account 
for 96 per cent of the total energy density of the universe 
(Spergel et al. 2003). An important prediction of CDM mod- 
els is the hierarchical formation of gravitationally bound 
structures. The smallest dark-matter haloes are expected 
to form first, whereas larger structures, ranging from spi- 
ral galaxies to galaxy clusters, are formed at a later stage 
through merging and accretion of smaller dark-matter haloes 
(bottom- up scenario). Numerical simulations of structure 
formation in CDM cosmologies have successfully reproduced 
the mass function and radial distribution of galaxies on the 
scales of galaxy clusters. On smaller scales, however, sim- 
ulations predict significantly more dark-matter haloes than 
being observed (Klypin et al. 1999; Moore et al. 1999). This 
discrepancy has been named the 'missing satellites' problem. 

To overcome this problem. Blitz et al. (1999) suggested 
that high-velocity clouds (HVCs) might be the gaseous coun- 
terparts of the 'missing' dark-matter haloes around the 



Milky Way. HVCs are gas clouds observed all over the 
sky in the 21-cm line emission of neutral atomic hydrogen. 
They were discovered with the Dwingeloo radio telescope by 
MuUer et al. (1963), and they are characterised by high ra- 
dial velocities of typically |wlsr| 100 kms~^ (see Wakker 
1991 for details). If HVCs were the 'missing' satellites, they 
could not have experienced significant star formation during 
their evolution. Therefore, they would appear in the form of 
pure gas clouds without any noticeable stellar population. 
In addition, HVCs would be spread all over the Local Group 
with typical distances of hundreds of kpc and high H I masses 
of about lO'^ Mq. 

At the same time, the expected large distances of HVCs 
from the Milky Way would result in fairly decent angular 
diameters of the clouds. Therefore, Braun & Burton (1999) 
defined a sub-sample of compact and isolated HVCs (the so- 
called CHVCs) which are characterised by angular sizes of 
less than 2° FWHM as well as isolation and separation from 
neighbouring Hi emission. The overall kinematics of the 
CHVC population is consistent with a distribution through- 
out the Local Group (de Heij et al. 2002a), making them 
promising candidates for the 'missing' dark-matter satellites 
around the Milky Way and the Andromeda Galaxy. Putman 
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Table 1. Observational parameters of the northern and southern 
part of our Effelsberg Hi blind survey of M31. The final baseline 
RMS is given for a system temperature of 25 K at the given 
velocity resolution. The specified sensitivity is the baseline RMS 
times the spectral bin width at the original velocity resolution, 
converted to H I column density. Sensitivities for the WNM were 
calculated for a brightness temperature of Tb > 3 cr, a velocity 
resolution of 20 kms~^, and a line width of 25 kms~^ FWHM. 



parameter 


south 


north 


unit 


autocorrclator 


old 


new 




bandwidth 


6.3 


10 


MHz 


polarisations 


2 


2 




channels per polarisation 


512 


4096 




velocity resolution 


2.6 


0.5 


kms~^ 


frequency switching 


normal 


in-band 




total integration time 


180 


180 


s 


final baseline RMS 


45 


60 


mK 


sensitivity 


21 


5.5 


10"' cm-2 


WNM sensitivity 


2.2 


1.3 


10^* cm-2 


WNM mass sensitivity 


8 


5 


10* Mq 



et al. (2002) extended the CHVC catalogue into the south- 
ern hemisphere, using the H I Parkes All-Sky Survey (Barnes 
ct al. 2001). The data for the northern and southern hemi- 
spheres were later combined by de Heij et al. (2002a) into 
an all-sky catalogue of 216 CHVCs. 

Several arguments have been raised against the idea of 
HVCs and CHVCs being the 'missing' dark-matter haloes 
predicted by CDM cosmologies. First of all, attempts to 
identify a population of HVCs in other nearby galaxy groups 
(Zwaan 2001; Braun & Burton 2001; Pisano et al. 2004) 
have failed, resulting in upjicr distance limits for HVCs and 
CHVCs from the Milky Way of the order of 150 kpc. Addi- 
tional evidence for HVCs being nearby at distances of the or- 
der of only 10 kpc has been provided by the detection of Ha 
emission from both HVCs and CHVCs (e.g., Kutyrev 1986; 
Weiner et al. 2001; Tufte et al. 2002; Putman et al. 2003) 
and from the determination of distance brackets for several 
HVC complexes (e.g., Danly et al. 1993; Wakker et al. 1996, 
2007,?; Thorn et al. 2006, 2008). Small distances from the 
Milky Way are also consistent with the head-tail structures 
found in numerous HVCs and CHVCs (e.g., Briins et al. 
2000, 2001; Westmeier ct al. 2005) which are thought to re- 
sult from ram-pressure interaction of the clouds with the 
ambient gas of the Galactic corona (Quilis & Moore 2001; 
Konz et al. 2002). 

The problem of determining the spatial distribution of 
HVCs can ultimately be solved by searching for the expected 
HVC jjopulation around the nearest large spiral galaxy, the 
Andromeda Galaxy (M31). First, the distance of M31 is well 
known so that important physical parameters of HVCs, such 
as their H I mass or their diameter, can directly be deter- 
mined from the observations. Second, we will look at the 
HVC population of M31 from the outside, allowing us to 
determine the (projected) radial distribution of HVCs. Fi- 
nally, M31 is relatively close to the Milky Way. Therefore, 
the sensitivity and angular resolution of large single-dish 
radio telescopes will easily allow us to detect HVCs of the 
mass and size of the large HVC complexes observed near the 
Milky Way. The first comprehensive search for HVCs around 



M31 was carried out by Thilker et al. (2004) with the 100-m 
Green Bank Telescope. They mapped an area of 7° x 7° in 
the 21-cm line of Hi with high sensitivity and discovered a 
population of about 20 HVCs out to the edge of their map at 
about 50 kpc projected distance from M31. Their discovery 
marked the first detection of an extensive HVC population 
around a galaxy other than our own. Several of these HVCs 
were studied in detail with the Westerbork Synthesis Radio 
Telescope (WSRT) by Westmeier et al. (2005). 

Thus, the radial extent of the population of HVCs and 
CHVCs around galaxies like the Milky Way and M31 is con- 
fined by two limits. An upper limit of about 150 kpc can 
be derived from the non-detection of HVCs in other galaxy 
groups, and a lower limit of about 50 kpc is marked by the 
edge of the H I survey of Thilker et al. (2004) out to which 
HVCs were found to be present. Therefore, we decided to 
carry out a complementary Hi survey with the 100-m ra- 
dio telescope at Effelsberg to search for HVCs and CHVCs 
around M31 out to much larger projected distances in excess 
of 100 kpc. This would allow us to trace the distribution of 
the HVC population of M31 over its entire radial extent and 
compare our results with the predictions of CDM structure 
formation scenarios. 

Our paper is organized as follows. In Sect. 2 we explain 
the technical aspects of our observations. In Sect. 3 our data 
reduction, calibration, and analysis strategy is described. In 
Sect. 4 we discuss the various completeness issues of our 
survey and their implications for our results. In Sect. 5 we 
describe the derived observational and physical properties of 
the high- velocity clouds detected in our survey. In Sect. 6 we 
discuss the results of our comparison of the observational pa- 
rameters of the HVCs near M31 with different CDM-based 
models and with the distribution of satellite galaxies around 
M31. Sect. 7 discusses the evidence for different hypotheses 
on the origin of HVCs. Finally, Sect. 8 summarises our re- 
sults and conclusions. 



2 OBSERVATIONS 

The observations for our H I blind survey of M31 were carried 
out between July 2003 and August 2004 with the 100-m 
radio telescope at Effelsberg. The observational parameters 
of the survey are summarised in Table 1. As a compromise 
between spatial coverage and observing time we chose the 
trapezoidal mapping area outlined in Fig. 1. It extends out 
to an angular distance of about 10° from M31 in the south- 
eastern direction, corresponding to a projected distance of 
about 140 kpc which is about two thirds of the projected 
distance towards M33. In this direction we expect the least 
confusion with foreground H I emission from the Milky Way. 
In the north-western direction our map extends out to about 
5°, corresponding to 70 kpc in projection. 

The map was split into individual rows of constant 
declination with 35 pointings each. The individual point- 
ings were separated in right ascension by 760 arcsec w 
HPBW X V2. Adjacent rows were separated in declination 
by 380 arcsec « HPBW/\/2 and shifted in right Eiscension 
by the same amount, resulting in an orientation angle of the 
entire map of 45° with respect to the equatorial coordinate 
system. As a consequence, the final map is beam-by-beam 
sampled and oriented almost perpendicular to the major 
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axis of the disc of M31. Aligning the map along the minor 
axis of M31 has the advantage of avoiding confusion with 
disc emission along the major axis as far as possible. 

The southern part of the map was observed with the old 
1024-channel autocorrelator, whereas for the northern part 
we used the new 8192-channel autocorrelator (AK90). The 
reason for changing from the old to the new autocorrelator 
in early 2004 was that in August 2003 the quality and sta- 
bility of the spectral baselines suddenly degraded. Several 
maxima and minima occurred over the bandpass, requir- 
ing a higher-order polynomial to approximate the shape of 
the baseline. The more severe problem was that the posi- 
tions and amplitudes of the maxima and minima changed 
significantly over short time scales of the order of one hour, 
making baseline correction difficult. As the quality of the 
spectra was not sufficient for the objectives of our M31 sur- 
vey we decided to carry out test observations with the new 
AK90. The test observations demonstrated that the AK90 
provided significantly better baseline quality and stability 
than the old 1024-channel autocorrelator. Therefore, we de- 
cided to continue the survey in the northern direction using 
the AK90 instead of the old autocorrelator. Another advan- 
tage of the AK90 is its larger bandwidth of 10 MHz. This 
allowed us to use the in-band frequency switching method 
which spends 100 per cent of the integration time on source. 
As a result the baseline RMS in the northern part of the 
map is by a factor of about \/2 lower than in the south- 
ern part. Throughout this paper wc follow the conservative 
approach of using the 1024-channel autocorrelator specifica- 
tions as the basis for all sensitivity-related parameters. 



3 DATA REDUCTION AND ANALYSIS 
3.1 Flux calibration 

The standard flux calibration source S7 was observed every 
6 hours and at the beginning and end of each observing run. 
The corresponding flux calibration factors for the two polar- 
isations were determined by the software nautocal (based 
on the results obtained by Kalberla et al. 1982), yielding 
statistical calibration errors of the order of only 1 per cent 
over an entire observing run of typically 10 to 15 hours. 



3.2 Stray radiation 

A potential problem could be that H i emission from the disc 
of M31 could have been detected through the near side lobes 
of the telescope and create the impression of extra-planar gas 
in proximity of the disc. The beam pattern of the Effelsberg 
telescope at A = 21 cm was studied in detail by Kalberla 
et al. (1980) through combined observations of Cygnus A 
and Cassiopeia A. The overall crosswise structure in their 
fig. 1 results from diffraction by the four support legs. Signif- 
icant side lobes can be found particularly in the north-south 
direction. At 1° angular distance from the pointing direction 
the sensitivity drops to about —30 dB. Typical brightness 
temperatures of 20 K, as observed along the H I ring of M31, 
would therefore be attenuated to about 20 mK which is far 
below the noise level of our data at the original spectral res- 
olution. The near side lobes in the east-west direction are 
not as extended as those in the north-south direction, but 
instead the sensitivity levels are somewhat higher near the 
pointing centre. At an angular distance of 15 arcmin the 
sensitivity drops to about —20 dB, resulting in an attenua- 
tion of a 20 K signal to about 200 mK. This is of the same 
order of magnitude as the expected signals of high-velocity 
clouds around M31. The near side lobes, however, will only 
affect regions of the map which are in projection close to 
the Hi disc of M31. All areas beyond a projected distance 
of the order of 10 kpc from the disc will not be influenced 
anymore by stray radiation from M31. In addition, stray ra- 
diation through the near side lobes should result in more 
extended emission, whereas most HVCs and in particular 
CHVCs are expected to be unresolved. Consequently, stray 
radiation will not affect our observations since most of the 
expected HVCs and CHVCs should be located at sufficiently 
large angular distances from the disc of M31 and will likely 
be unresolved by the HPBW of 9 arcmin. 



3.3 Spectral bsiseline correction 

To improve the quality of the spectral baseline we imple- 
mented a two-step baseline correction procedure. The first 
step was applied to each of the two polarisations separately. 
All 35 spectra in a single row of the map were averaged, and 
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Figure 2. Spatial completeness of the M31 survey as a function 
of projected distance, r, from the centre of M31. The dashed line 

shows the global completeness function, CQ{r), and the solid line 
the local completeness function, C'i^{r). The vertical dotted lines 
indicate the radii a, 6, c, and d at which the completeness function 
changes abruptly due to the special geometry of the observed field. 
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signal-to-noise ratio 

Figure 3. Fraction of detected spectral lines as a function of 
signal-to-noise ratio, derived from the investigation of 500 arti- 
ficial spectra as described in the text. About 50 per cent of all 
signals with a peak intensity of 3.5 tr were detected. Above 5 a we 
arc virtually complete with only a single of the 14 lines above 5 a 
being missed. 



a polynomial of usually S**" order over Av = 900 kms~^ was 
fitted to the averaged spectrum after manually setting line 
windows. This constant reference polynomial - representing 
the average baseline over a time period of about two hours 
- was then subtracted from each individual spectrum of the 
row. The resulting spectra were already quite smooth as the 
baseline shape at Effelsberg was fairly stable over a period 
of two hours. The second step was to average all spectra at 
a single position on the sky and fit a polynomial to the re- 
sulting spectrum. Because of the previous subtraction of the 
reference polynomial, a low-order polynomial WEis sufficient 
for the final averaged spectrum. In most cases, the baseline 
shape could be sufficiently described by a 4**^ or 5**^ order 
polynomial across a velocity range of 900 kms"^ Only m 
very few cases a lower (down to I''') or higher (up to 8**^) 
order was necessary. After cutting out a few strong radio 
frequency interference signals by hand, the resulting spectra 
were Hanning-smoothed to 4.1 kms~^ velocity resolution 
and combined into a FITS data cube for further analysis. 

3.4 Identification of sources 

The final survey consists of 3586 individual pointings for 
each of which an H I spectrum in the velocity range of vlsr ~ 
—750 . . . -1-150 kms~^ was available for further analysis. As 
we were interested in the detection of HVCs and CHVCs 
around M31, the survey was particularly optimised for the 
detection of compact sources with narrow spectral lines. lu 
the in-band frequency switching mode nearly 100 per cent 
of the integration time is spent on source, maximising the 
sensitivity for the detection of faint CHVCs. The spectral 
baselines in this special observing mode, however, are not 
particularly smooth and stable, restricting the detectability 
of faint, diffuse sources. 

To identify HVC/CHVC candidates of M31 all spectra 
were searched by eye for potential emission lines. The cri- 
terion for HVC emission was a separation of the signal in 
phase space from the disc emission of M31 (and the Milky 
Way), where phase space is considered the observable three- 



dimensional sub-space consisting of the two spatial coordi- 
nates in the plane of the sky and the radial velocity as the 
only available kinematic component. Thus, signals could ei- 
ther be located on the same line of sight as the disc emission 
of M31 but at different radial velocities, or they could be lo- 
cated outside the disc at any radial velocity. At the same 
time, we did not consider objects with v-lsr > —140 kms""' 
to avoid confusion with Galactic foreground emission. 

It is important to note that we have to change the defi- 
nition of 'high-velocity cloud' for our M31 observations. The 
HVCs around the Milky Way were characterised by their pe- 
culiar radial velocities as seen from a point of view within 
the Galactic Hi disc. In the case of M31, however, we are 
looking at the entire galaxy and its HVC population from 
the outside. Therefore, the term 'HVC refers to all circum- 
galactic gas clouds which are thought to be the equivalents of 
the circumgalactic gas clouds around the Milky Way known 
as 'high-velocity clouds'. This includes clouds whose radial 
velocities are compatible with the rotation curve of M31 but 
which are located well beyond the edge of the H I disc. 



4 COMPLETENESS OF THE DATA 
4.1 Spatial coverage of the map 

An important question arising from the approach of a vi- 
sual inspection of the survey is the question of detectability 
and completeness of HVC/CHVC candidates. A first limi- 
tation of the completeness of our M31 survey results from 
the spatial incompleteness of our map at larger radii. The 
geometric situation is illustrated in Fig. 1. There are five dif- 
ferent regions for which analytic expressions for the spatial 
completeness have to be derived separately. The five regions 
are separated from each other at radii of 2?5, 3?5, 7?9, and 
12?7 from the centre of M31, corresponding to projected dis- 
tances of 34 kpc, 48 kpc, 108 kpc, and 174 kpc if the distance 
of M31 is assumed to be 780 kpc (Stanek & Garnavich 1998). 
There are two ways of calculating the spatial completeness 
within these regions as a function of radial distance from the 
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Figure 4. Hi column density maps of our Effclsbcrg blind survey of M31. Top: Integrated Hi column density in the velocity range of 
''LSR = —620...— 24 kms^'^. The contours were drawn at 4, 8, 16, and 25 X lO'^" cm~-^. The faint, filamentary emission all over the 
map originates from Galactic Hi. Bottom: Integrated Hi column density in the velocity range of ?;lsr = —620 . . . —139 kms"'^, using a 
logarithmic intensity scale. The black contours correspond to 2 x 10^* and 1 x 10^^ cm~^. The white contours are those from the upper 
map. Several regions of extra-planar gas and high-velocity clouds can be seen all around M31. 



centre of M31. One way is to calculate the local filling factor 
at distance r from the centre of M31 which is given by the 
ratio of the length of the arc intersecting the map over the 
total arc length, 2nr, of the full circle. The other way is to 
calculate the global filling factor of the map within a certain 
radius, r, which is the ratio of the area being covered by the 
map within r over the total area, tt . 

The two spatial completeness functions, Ci^ir) and 
Co(r), are plotted in Fig. 2 as a function of projected dis- 
tance, r, from the centre of M31. Within r = 34 kpc the 
survey is complete. Beyond r = 174 kpc the local complete- 
ness drops to whereas the global completeness decreases 
with At intermediate radii both functions have a more 

complex behaviour. They drop abruptly at r = 34 kpc with 
the local completeness function, Ci^{r), decreasing much 
faster than the global completeness function, Coir). Lo- 
cally, 50 per cent completeness is reached at a distance of 
T'50 = 48 kpc already, whereas on the global scale the cor- 
responding radius is rso = 77 kpc. At a projected distance 
from M31 of 100 kpc the completeness is Cl = 12 per cent 
and Cg = 36 per cent. This means that out to 100 kpc our 



survey still covers more than a third of the total area, but 
the local completeness at this distance is already fairly low 
because the survey area is concentrated around the central 
region. 



4.2 Spatial sampling of the map 

Another reason for spatial incompleteness is related to the 
spatial sampling of the map. The individual pointings of 
the map are separated by the HPBW of the telescope of 
9 arcmin, so that the sensitivity is significantly decreased at 
positions which are located exactly in between four individ- 
ual pointings. Assuming a two-dimensional Gaussian sensi- 
tivity function for the main beam, the sensitivity of each 
beam at the position in between four neighbouring beams 
(i.e. the corner of each grid element) is 0.25. Averaging the 
four neighbouring beams results in a decrease of the RMS 
by a factor 2 and, thus, in a slightly higher effective sensitiv- 
ity of 0.5. As a consequence, the sensitivity for unresolved 
sources is in the range of 0.5 to 1, depending on the source 
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Figure 5. Left: Locations of the 17 individual HVCs identified in our Eff'clsbcrg survey (filled circles) superposed on an Hi column 
density map. Two HVCs studied with the WSRT (Westmeier et al. 2005) but not detected in our survey are plotted as solid open circles. 
The dashed open circles denote the positions of known satellite galaxies of M31. Right: Hi column density map in the velocity range 
of ^'LSR = —608 . . . —135 kms~^, using the same mask as for Fig. 4. The contours range from 18 to 22 dex in steps of 0.5 dex. The two 
dashed arrows indicate the lines for the position-velocity diagrams in Fig. 6. 



position relative to the grid. Extended sources are covered by 
several pointings so that their detectability is not reduced. 



4.3 Flux detection limit 

Another important question is down to which signal-to-noise 
ratio potential CHVCs can be detected. To answer this ques- 
tion, we simulated 500 spectra with a Gaussian RMS and a 
velocity resolution of 10.1 kms^^. Each spectrum was multi- 
plied with a sine function of random amplitude, wavelength, 
and phase shift to get arbitrary baseline shapes. With a 
probability of 15 per cent a Gaussian spectral line was added 
to each spectrum. Line widths were randomly chosen be- 
tween 15 and 35 kms~^, corresponding to what is expected 
for the warm neutral medium of HVCs. Peak intensities were 
also assigned randomly with signal-to-noise ratios between 
and 6. Our algorithm added 81 artificial spectral lines with 
the above-mentioned parameters to the 500 spectra. Posi- 
tions, intensities, and line widths of the 81 spectral lines 
were written to a data file for later comparison. 

Next, the entire data set of 500 artificial spectra was 
searched by eye to detect the hidden spectral lines. In to- 
tal, 34 of the 81 signals could be identified. The resulting 
completeness function is shown in Fig. 3. At a peak inten- 
sity level of 3.5 a the completeness is about 50 per cent. 
Above 5cr we are virtually complete with only a single of 
the 14 lines above 5 a being missed. Note that there arc 
also two detections at about 1 cr. In both cases strong noise 
peaks by chance have increased the signal above the detec- 
tion threshold. These results demonstrate that, when using 
the classical 3 <t criterion for the detectability of spectral sig- 
nals, we have to consider that source counts below the 5 a 
level suffer from incompleteness. On the other hand, detec- 
tions are possible even at intensities below 3 a although the 
detection rate decreases significantly. 



4.4 Confusion with Galactic emission 

In optical astronomy the investigation of cxtragalactic ob- 
jects is affected by extinction due to dust in the interstellar 
medium of the Galaxy. The strongest extinction is observed 
at low Galactic latitudes along the so-called 'zone of avoid- 
ance'. The situation in Hi spectroscopy is slightly different. 
Emission lines of cxtragalactic objects can be superposed 
on those of the Milky Way, making an identification and 
analysis difficult or even impossible. Galactic Hi emission 
is found not only near the Galactic plane but all across 
the sky. However, only the radial velocity range around 
vlsr = kms^^ is affected, whereas sources with radial 
velocities of [vlsrI ^ 100 kms~^ are usually not subject to 
confusion with Galactic emission lines. In other words, the 
'zone of avoidance' in H I spectroscopy covers the entire sky 
but only a limited radial velocity range for each line of sight. 

The systemic radial velocity of M31 is «lsr ~ 
—300 kms~^. On the approaching side of M31 we observe 
radial velocities down to about —600 kms~^, implying val- 
ues up to about krris^^ on the receding side. The latter is 
well within the velocity range of Galactic H l emission in this 
area which reaches out to about —140 kms~^. Hence, a mar 
jor part of the receding side of M31 is blended with Galactic 
emission. Blending effects could also affect the detection of 
HVCs around M31. As the peak intensities of HVCs are ex- 
pected to be small compared to the intensities of Galactic 
emission lines, HVCs superposed on Galactic emission are 
virtually undetectable. In our survey this should be the case 
for all potential HVCs with wlsr > —140 kms^^, resulting 
in incompleteness with respect to certain radial velocities. 

5 THE DETECTED POPULATION OF HVCS 

In total we identified 17 discrete HVCs in our Effelsberg 
survey. Each HVC was later confirmed by follow-up observa^ 
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Figure 6. The two panels show position-velocity diagrams along 
the lines indicated by the two dashed arrows in the column density 
map in the right panel of Fig. 5. Several HVCs and regions of 
extra-planar gas can be seen around M31. The contour levels 

arc 50, 100, and 200 mK. Note that Davics' Cloud is completely 
isolated in phase space although in the column density map there 
appears to be a connection with the Hi disc of M31. 



tions with the Effelsberg telescope. An overview of all clouds 
and their observational parameters is given in Table 2 and 
Fig. 4 and 5. All HVCs are located within a projected dis- 
tance of about 50 kpc from M31. Some of them are isolated 
in position and velocity, whereas others appear to be asso- 
ciated with the diffuse extrarplanar gas near the edge of the 
Hi disc of M31. A few clouds, such as number 14 in our 
catalogue, arc seen in projection against the Hi disc of M31 
and can only be discriminated by their different radial ve- 
locities. This situation is illustrated by the position-velocity 
diagrams in Fig. 6. Hi spectra of all clouds are shown in 
Fig. 7. 

Hi column density maps of the entire survey area 
are displayed in Fig. 4. The upper map shows the en- 
tire flux integrated over a velocity range of vlsr = 
-620 . . . -24 kms"^ In the direction of M31 the Hi emis- 
sion of the Milky Way disc reaches out to velocities of 
vlsr ^ —140 kms~^, resulting in a significant overlap with 
emission from the north-eastern part of M31. Filamentary 
emission from the Galactic disc can be seen all over the map 
with typical column densities of a few times lO'^" cm~^. Two 
conspicuous features of M31, the famous Hi ring and the 
warp of the H l disc, are immediately visible in the map. The 
ring was first clearly detected by Roberts (1966) and has a 
radius of approximately 10 kpc. The central depression in 
the Hi disc shows column densities of about 30 per cent of 
the peak values found along the ring (Guibert 1974). The 
warping of M31 can be seen towards the far ends of the disc, 
in particular in the north-eastern part, where the so-called 



'perturbed outer arm' (Whitehurst et al. 1978; Sawa & Sofue 
1982) is significantly bent in the northern direction. 

The lower map in Fig. 4 is a high-contrast logarithmic 
map of H I column densities of 10^* . . . 10^^ cm~^ in the ve- 
locity range of ulsr ~ —620 . . . —139 kms~^. Here, the ve- 
locity range of Galactic H i emission has been excluded with 
the restriction that some of the emission from the north- 
eastern end of the disc of M31 is also missing. In addition, 
we created a mask data cube which was smoothed in posi- 
tion and velocity. Only those mask elements with intensities 
in excess of 17.5 mK were selected in the original data cube 
to contribute to the moment map, resulting in a high- 
quality image showing Hi column densities down to a few 
times 10^** cm"^ 

Apart from regular disc emission, several regions of 
extra-planar gas and isolated Hi clouds are visible. Extra- 
planar gas emission is particularly prominent along the 
south-eastern edge of the disc, but extended, diffuse H I 
emission with column densities of Nu_\ < 10^® cm~^ is also 
present along the north-western edge. The most outstanding 
isolated H i cloud in the map is Davies' Cloud located about 
1?5 north-west of the centre of M31. The cloud was discov- 
ered by Davies (1975) and studied in detail by de Heij et al. 
(2002b) with the WSRT. A few more compact and isolated 
clouds can be seen all around M31. 

Most of the clouds seen in Fig. 4 had already been dis- 
covered by Thilker et al. (2004) in their H I survey of M31 
with the Green Bank Telescope. A comparison between their 
fig. 2 and the lower map in Fig. 4 reveals a good correla- 
tion between the clouds seen in both surveys with the ex- 
ception of a few faint and extended clouds detected in the 
GET data but not seen in the Effelsberg map. Those clouds 
without counterparts in the Effelsberg survey are typically 
very diffuse with H I column densities of the order of only 
10^^ cm~^. They were only detected in the GET data after 
decreasing the velocity resolution significantly to 72 kms^^. 
With the frequency switching method used for our Effelsberg 
observations such diffuse objects can hardly be detected. In- 
stead, our survey was optimised for the detection of compact 
sources with narrow spectral lines. 

Within the 3 a mass detection limit of our survey of 
Mhi ~ 8x 10* Mq we did not find any extended population 
of hundreds of CHVCs as observed around the Milky Way. In 
particular, we did not detect any HVCs or CHVCs beyond a 
projected distance of about 50 kpc from M31, although our 
survey reaches out to significantly larger projected distances. 
Apparently, the detected HVCs are concentrated near M31 
and less numerous than expected from the population of 
hundreds of HVCs/CHVCs observed around our Milky Way. 



5.1 Parameters of the detected HVCs 

5.1.1 Velocities 

The radial velocities of the HVCs are in the range of 
«LSR ~ —520 . . . —150 kms~^. The distribution of radial ve- 
locities is not perfectly symmetric with respect to the sys- 
temic velocity of M31 of vlsr « -300 kms"^ (Fig. 8). First 
of all, confusion with the H I emission of the Milky Way 
does not allow us to detect any HVCs at velocities above 
— 140 kms~^, explaining the upper limit in the observed ve- 
locity range. Therefore, one or two more HVCs beyond this 
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Table 2. List of the 17 discrete HVCs around M31 detected in our H i survey with the Effelsberg telescope. The columns denote the 
catalogue number of the cloud, its name according to Westmeier et al. (2005), right ascension, a, and declination, S, the projected 
distance, ^proj, of the column density maximum from the centre of M31, the LSR radial velocity, i^lsRi the FWHM of the spectral line, 
Av, the peak Hi column density, AThi, and the Hi mass, Mhi- Radial velocities are accurate within about ±5 kms~^, and the statistical 
uncertainties of line widths, column densities, and H i masses are of the order of 10 per cent. 



# 


Name 


a 


5 


''proj 


■vlsr 




A^HI 


Mm 






(J2000) 


(J2000) 


(kpc) 


(kms-i) 


(kms^"'") 


(10^'' cm-2) 


(lOS Mo) 


1 


M31 HVC 1 


00'^43™00" 


38° 53' 


32.5 


-501 


17 


1.1 


7.8 


2 


M31 HVC 2 


00^41^06" 


39° 16' 


27.6 


-509 


34 


1.7 


6.3 


3 


M31 HVC 7 


00^410106= 


39° 28' 


24.9 


-463 


47 


1.2 


3.9 


4 


M31 HVC 5 


Ooh43moos 


39°31' 


23.9 


-507 


34 


0.9 


6.3 


5 


M31 HVC 12 


00''43™17^ 


39°35' 


23.0 


-409 


58 


1.1 


4.3 


6 


M31 HVC 10 


00''41™05= 


39°41' 


22.0 


-426 


26 


0.9 


4.7 


7 




00'^33™18= 


40°32' 


26.2 


-518 


39 


0.7 


3.1 


8 


Davies' Cloud 


00^37^25= 


42° 26' 


20.7 


-442 


21 


3.9 


130 


9 




00h47™43= 


40° 19' 


18.3 


-324 


71 


1.6 


5.9 


10 




00^47^12= 


40°51' 


12.9 


-282 


67 


1.4 


5.0 


11 




00'^43™53= 


42° 07' 


11.9 


-266 


16 


0.4 


1.2 


12 




00''45™35= 


42°00' 


12.3 


-223 


13 


0.3 


1.0 


13 


M31 HVC 15 


00''43™55^ 


44° 13' 


40.2 


-266 


30 


0.9 


3.5 


14 


M31 HVC 13 


00''39™56^ 


41° 16' 


7.1 


-207 


14 


0.5 


6.4 


15 




00''45™02= 


42°45' 


21.0 


-206 


21 


0.5 


1.7 


16 




00^39^09= 


45° 04' 


52.5 


-165 


17 


0.6 


2.2 


17 




00'^ 32™ 42= 


44° 20' 


48.7 


-149 


11 


0.9 


8.2 



velocity limit could be missing in our sample due to confu- 
sion with Galactic omission. On the other hand, several of 
the detected clouds are arranged in a complex of HVCs near 
the south-eastern edge of the disc of M31 with very similar 
radial velocities in the range of about —500 . . . —400 kms~^. 

The possibility of confusion with Galactic Hi emission 
immediately raises the question whether HVCs near Galac- 
tic velocities might actually be associated with the Milky 
Way rather than M31. The detection of very compact and 
faint gas clouds and HVCs in the halo of the Milky Way 
(e.g., Lockman 2002; Briins & Westmeier 2004; Hoffman 
et al. 2004; Richter et al. 2005) has shown that the criteria 
of small angular size and low H I column density axe not suf- 
ficient to discriminate between HVCs near the Galaxy and 
M31. Given the lack of distance information, it is in principle 
impossible to associate each individual cloud with M31. The 
observational parameters of the detected HVCs, however, al- 
low us to definitely associate the population as a whole with 
the Andromeda Galaxy. First of all, the clouds are clearly 
concentrated around M31 with projected distances of less 
than about 50 kpc. If they were associated with the Milky 
Way we would expect them to be more equally spread across 
the survey area. The probability to have them distributed 
in a 50 kpc circle around M31 by chance is of the order 
of only 10~® in the case of a random distribution across 
the sky. Second, the observed radial velocities are nearly 
uniformly distributed around the systemic velocity of M31 
of —300 kms~^, although we arc probably missing one or 
two HVCs with velocities above —140 kms^^. Furthermore, 
those HVCs with velocities of about —500 kms~^ would 
have significantly higher negative velocities than almost any 
other HVC near the Milky Way. Therefore, both the spatial 
and kinematic parameters of the detected HVCs strongly 
suggest that the population as a whole is associated with 
M31, although we cannot make this statement for individ- 
ual clouds. 



5.1.2 Line widths 

The observed line widths of the HVCs cover a very large 
range from about 10 to 70 kms"^ FWHM (Fig. 8). Line 
widths of less than 20 kms^^ indicate the possible pres- 
ence of a cold neutral gas component or a multi-phase 
cold/ warm neutral medium which is probably unresolved 
by our 9 arcmin HPBW. A multi-phase medium is com- 
monly found in Galactic HVCs and CHVCs (see, e.g., Cram 
& Giovanelh 1976; Briins et al. 2000, 2001; Westmeier et al. 
2005). It can also be reproduced in hydrostatic simulations 
of HVCs (Sternberg et al. 2002) so that our evidence for the 
presence of a multi-phase medium in HVCs around M31 is 
not surprising. Difficult to explain, however, is the detec- 
tion of extremely broad lines in some of the clouds. In a few 
cases line widths as large as 60 or 70 kms~^ FWHM are 
observed which can no longer be explained by thermal line 
broadening alone as the required temperatures would exceed 
10^ K. The most likely explanation for large line widths in 
excess of 30 or 40 kms"'^ FWHM is a high degree of internal 
kinematics within the area covered by the 9 arcmin HPBW 
of the Effelsberg telescope (about 2 kpc in projected size 
at the distance of M31). On the one hand, several clouds 
or clumps with slightly different radial velocities could be 
aligned along the line of sight. On the other hand, individ- 
ual clouds could have regular or turbulent internal motions 
such as global rotation or gas flows. 



5.1.3 Column densities and masses 

The detected H i column densities are typically of the order 
of 10^^ cm~^. This very low value indicates that most of the 
HVCs are probably not resolved by the Effelsberg telescope 
so that the column densities listed in Table 2 have to be con- 
sidered lower limits averaged over the 9 arcmin HPBW of the 
telescope. In fact, our high- resolution synthesis observations 
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Figure 7. Effelsberg H I spectra of the HVCs around M31 according to Table 2. The position of the spectral line belonging to the HVC 

is marked with a vertical line in each spectrum. The strong emission at vlsr ~ kms~^ is usually due to Hi gas in the Galactic disc. 
In a few spectra, however, emission from the M31 disc is also present at different velocities. 
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Figure 8. Histograms of radial velocities (top) and line widths 
(bottom) of the 17 individual HVCs identified in our survey. 



of some of the HVCs with the WSRT (Westmeier et al. 2005) 
revealed significantly higher peak column densities compa- 
rable to what is observed in Galactic HVCs. Assuming a 
distance of 780 kpc (Stanek & Garnavich 1998), the result- 
ing Hi masses of the HVCs are typically in the range of a 
few times 10'^ Mq with a mean value of 4.5 x 10'^ Mq. This 
excludes Davies' Cloud which - if at the distance of M31 - 
would have a significantly higher Hi mass of 1.3 x lO'^ M©. 
In fact, Davies' Cloud would contribute about 65 per cent of 
the total mass of all detected HVCs of 2 x lO'' M©. The total 
Hi mass of all individual HVCs and regions of extra-planar 
gas around M31 together is of the order of 5 x 10^ Mq. The 
true value will be slightly higher as some of the gas may be 
missing due to blending with the H I emission of the Milky 
Way and the disc of M31 itself. 

It is interesting at this point to compare the H I masses 
of the HVCs around M31 with those of Galactic HVC com- 
plexes. For several complexes, namely Complex M (Danly 
et al. 1993; Ryans et al. 1997), Complex A (Wakker et al. 
1996; van Woerden et al. 1999), Complex WB (Thom et al. 
2006), Complex WE (Wakker 2001), Complex C (Wakker 
et al. 2007; Thom et al. 2008), and the Cohen Stream and 
Complex GCP (Wakker et al. 2007), distance brackets have 
been published so far. These complexes are typically located 
a few kpc above the Galactic plane, implying H i masses in 
the range of log(MHi/M0) w 4.8... 6. 7. This mass range 
is consistent with the Hi masses observed for the HVCs 



around M31, supporting our conclusion that they represent 
the analogues of the large HVC complexes seen around the 
Milky Way. Only the Magellanic Stream and the Leading 
Arm have significantly higher masses. Due to the lack of dis- 
tance information, their masses can only be estimated from 
numerical simulations in connection with Hi observations. 
Connors et al. (2006) derived H I masses of 2.4 x 10® M© and 
7.3 X lO'^ Mq for the Magellanic Stream and the Leading 
Arm, respectively, which are by orders of magnitude higher 
than those of the HVCs around M31, including Davies' 
Cloud. The Magellanic Clouds, however, are a special case 
with no equivalent in the M31 system. 



5.1.4 Dynamical masses and external pressure 

To get an estimate of the mass fraction of neutral hydrogen 
in the M31 HVCs we can compare the observed H I masses of 
the HVCs with the FWHM of the spectral lines which, in the 
ideal case, should provide an estimate of the virial masses 
of the clouds. In Fig. 9 we have plotted the Hi mass, Mhi, 
versus the FWHM, Av, of all HVCs detected in our sample. 
The error in line width is assumed to be 10 per cent, whereas 
the error bars for the mass indicate our 3 a mass sensitivity 
of 8 X 10* Mq . If their H I masses were proportional to their 
virial masses, all clouds would be scattered around a straight 
Hue running from the lower left to the upper right part of the 
diagram. In its most simple form the virial equation reads 



^2 8 In 2 GMvir GMvir 
Av = Rs 

5 Rvir Rvir 



(1) 



assuming an isothermal sphere with constant mass density 
and the same mass for all particles forming the cloud. We 



can now replace the virial mass 
mass by introducing a factor a = 



Mvir, in Eq. 1 by the Hi 



Av'' 



Mvir /Mhi such that 
GMhi 

RyiT 



(2) 



Lines of constant virial-to-Hl mass ratio, a, according to 
Eq. 2 are also plotted in Fig. 9, assuming a constant virial 
radius of -Rvir = 0.5 kpc for each cloud. This value for R^ir 
is suggested by our follow-up observations of several of the 
HVCs with the WSRT (Westmeier et al. 2005). Obviously, 
most of the HVCs near M31 are spread around a w 100, 
suggesting that they could be gravitationally dominated by 
other mass components such as ionised gas or dark mat- 
ter. However, a few clouds show very high mass ratios of 
a « 1000. In most of these cases we observe exceptionally 
broad H I lines with line widths exceeding 40 km s~^ FWHM, 
suggesting that the clouds cannot be considered entities in 
dynamical equilibrium so that the observed line widths prob- 
ably do not provide any meaningful estimate of the asso- 
ciated total mass. On the other hand, Davies' Cloud lies 
slightly off the other data points, resulting from its fairly 
narrow spectral lines combined with a very high H I mass. 
None the less, its true mass ratio at the distance of M31 
would be a « 26 because its physical diameter of about 
7 kpc would be much higher than the value of 0.5 kpc as- 
sumed in our calculations. It is important to keep in mind, 
however, that for all these calculations we had to assume 
that the clouds are at the same distance as M31. 

Another possibility to stabilise the clouds in addition 
to their own mass is through the pressure of their ambient 
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Figure 9. Hi mass, M"hi, plotted against the FWHM of the 
spectral line, At!, for all HVCs detected in our Effelsbcrg survey. 
The dotted Unes indicate different values of a = M^\^ / M^n under 
the assumption of a constant virial radius of Rvir = 0.5 kpc (see 
Eq. 2). Hence, clouds with a small M^^^/Myh ratio are located in 
the upper left region of the diagram, whereas clouds with a large 
Mvir/Mni ratio can be found at the lower right. 



medium. Assuming an isothermal sphere, the virial equation 
for a cloud with external pressure support reads 

3kTMvir 3GM^i, 



(3) 



/i 5-Rvir 

where P is the external pressure, T is the g£is temperature, 
and fj, denotes the mean mass per particle (Spitzer 1978). 
Let us assume the virial mass to be equal to the Hi mass 
and of the order of Mvir = 10'' M0. Then we can use Eq. 3 
to determine the external pressure required to stabilise the 
cloud without any additional mass component. Assuming a 
virial radius of Rvir = 0.5 kpc, a temperature of the warm 
neutral medium of T = 10* K, and a moan mass equal to 
the mass of a hydrogen atom of /i = 1.674 x 10^^^ kg, we 
obtain an external pressure of 

P/k = 84Kcm"^. (4) 

Therefore, an external pressure of the order of 100 Kcm~^ 
is sufficient to stabilise a typical HVC around M31 in addi- 
tion to its own gravitational potential. In this case uo ad- 
ditional mass components, such as molecular gas or dark 
matter, would be necessary. (Note that our calculation does 
not even account for the expected amount of primordial 
helium.) The required pressure is of the order of what is 
expected for the ionised coronal gas around massive spiral 
galaxies with a temperature of about 10® K and a typical 
density of 10"'' cm~^ (e-g-, Baldwin 1954; Sembach et al. 
2003; Rasmusson ct al. 2003). This result demonstrates that 
the HVCs observed around M31 can be in pressure equilib- 
rium with the circumgalactic hot corona without the need 
for any additional dark matter component for stabilisation. 

5.2 Discussion of individual objects 

5.2.1 Davies' Cloud 

As mentioned before, Davies' Cloud would be the by far 
most massive individual HVC, if it were at the same dis- 



tance as M31. Discovered by Davies (1975), it was studied 
in detail with the WSRT by de Heij et al. (2002b). It is lo- 
cated only about 1?5 north-west of M31, corresponding to a 
projected distance of only 20 kpc. With an angular diame- 
ter of about 30 arcmin FWHM, corresponding to a physical 
size of roughly 7 kpc at the distance of M31, Davies' Cloud 
is also the most extended HVC near M31. Its appearance is 
not perfectly symmetric, but the column density maximum 
is slightly shifted in the direction of M31 (see Fig. 10 for a 
column density map and example spectrum). The resulting 
head-tail morphology could be an indication for a distortion 
of Davies' Cloud by the ram pressure of its ambient medium. 
In the column density map in Fig. 5 Davies' Cloud appears 
to be connected with the disc emission of M31. However, 
there is no such connection in phase space because the gas 
in the putative 'bridge' has completely different radial veloc- 
ities and does not provide any connection between Davies' 
Cloud and M31. This isolation of Davies' Cloud in phase 
space is obvious from the position-velocity diagram in the 
right panel of Fig. 6. 

Whether Davies' Cloud is indeed associated with M31 
or whether it is a foreground object in the vicinity of the 
Milky Way has been discussed in some detail by dc Heij 
ct al. (2002b). In their high-resolution synthesis observa- 
tions with the WSRT they resolve the cold gas component 
of Davies' Cloud into an extended arc of compact clumps 
stretching mainly across the eastern edge of the cloud. This 
configuration, in connection with the asymmetry of the dif- 
fuse neutral gas component, suggests that Davies' Cloud is 
interacting with an ambient low-density medium, such as 
the gaseous halo or corona of M31. 

However, another possibility, originally discussed by 
Davies (1975), is an association of Davies' Cloud with the 
Magellanic Stream. Braun & Thilker (2004) detected a faint 
northern extension of the Magellanic Stream in their WSRT 
total-power survey of a large area in the region around 
M31 and M33. In fact, Davies' Cloud lies within only a 
few degrees of parts of the stream. Although typical col- 
umn densities of the stream in this region are much lower 
than that of Davies' Cloud, the radial velocity of Davies' 
Cloud (ulsr = —442 kms~^) is very similar to the ve- 
locities of the gas associated with the Magellanic Stream 
(wlsr — —450 kms~^). Therefore, a connection of Davies' 
Cloud with the Magellanic Stream is plausible, and Davies' 
Cloud may be a foreground object in the vicinity of the 
Milky Way. Consequently, Davies' Cloud may be less mas- 
sive than previously assumed. Unfortunately, the distance of 
the far end of the Magellanic Stream is not well constrained, 
but if we consider a distance of the order 50 kpc for Davies' 
Cloud, ^ its Hi mass would only be about 5 x 10* Mq. 

In fact, there is more evidence for Davies' Cloud being 
part of the Magellanic Stream. At the distance of M31 the 
Hi mass of Davies' Cloud would be 1.3 x lO'^ M0, which is 
comparable to the largest HVC complexes near the Milky 
Way. However, while most of the HVC complexes near the 
Milky Way have a very elongated and filamentary structure, 
Davies' Cloud appeaxs more or less spherical. This morpho- 



^ Recent Arecibo observations by Stanimirovic et al. (2008) sug- 
gest a distance of the order of 70 kpc for the far northern end of 



the Magellanic Stream. 
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Figure 10. Hi column density map of Davics' Cloud. The outermost contour corresponds to 2 X 10 cm~ , and the remaining contour 
lines are drawn in steps of 5 X 10^** cm^"^, starting at 5 X lO'^* cm^'^. On the right hand side the Hi spectrum at the position of the 
column density maximum is shown. With an Hi mass of 1.3 x lO'^ Mq Davies' Cloud would be the by far most massive HVC near M31. 



logical difference is surprising. It is very unlikely that the 
cloud is significantly elongated along the line of sight as in 
this case we would expect to sec very complex and broad 
spectral line profiles caused by a complex spatial and kine- 
matic structure along the filament. Instead, the observed Hi 
lines are rather narrow with Aw « 20 kms"^ (FWHM) and 
of Gaussian shape (Fig. 10), confirming our impression of 
an overall spherical morphology of Davies' Cloud. 

Instead, the morphology of Davies' Cloud is reminiscent 
of compact HVCs near the Milky Way and in particular the 
numerous compact clouds in the vicinity of the Magellanic 
Stream and the Leading Arm. Ben Bekhti et al. (2006) stud- 
ied two of these clouds near the Leading Arm with Parkes 
and the ATCA. They detected numerous small clumps in 
the heads of the clouds similar to those found in Davies' 
Cloud by de Heij ot al. (2002b). In fact, the similarity be- 
tween fig. 4 of Ben Bekhti et al. (2006) and fig. 2 of de Heij 
et al. (2002b) is compelling. The H I lines observed in Davies' 
Cloud have significantly laxger line widths, but this can eas- 
ily be explained by a superposition of clumps along the line 
of sight or by slightly different physical conditions, all the 
more so as relatively broad H I lines are commonly observed 
in the cores of compact HVCs (e.g., Braun & Burton 2000). 

Therefore, we have to conclude that Davies' Cloud is 
most likely associated with the Magellanic Stream and not 
located at the distance of M31. This is suggested by its ra- 
dial velocity, its physical properties, and its morphological 
structure. Irrespective of its origin, we include Davies' Cloud 
in our analysis of M31 HVCs for the sake of completeness, 
albeit emphasising its special status where indicated. 

5.2.2 Extra-planar gas south-east of M31 

One of the most outstanding features detected in our survey 
is the region of extra-planar gas near the south-eastern edge 
of the Hi disc of M31. To study this area with better sen- 
sitivity we re-observed it with a longer integration time of 
5 minutes per position in in-band frequency-switching mode 
with the Effelsberg telescope, using the new autocorrclator 
AK90. The parameters are otherwise similar to those of our 
blind survey. The resulting H I column density map of the 
entire region is shown in the left panel of Fig. 11. The to- 



tal Hi mass of the extras-planar gas in this area is about 
2 X 10^ M0, thus even exceeding the Hi mass of Davies' 
Cloud if it were at the distance of M31. 

From the map it is obvious that the extra-planar gas is 
not totally diffuse but instead forming more or less compact 
concentrations and clumps. The clumpy structure is spar 
tially not well resolved by our 9 arcmin HPBW. However, 
the individual clumps arc characterised by different radial 
velocities which are kinematically well resolved by our spec- 
tral resolution. Two example spectra are presented in the 
right panel of Fig. 11. In many cases the spectral profiles are 
rather complex and vary significantly between neighbouring 
pointings of the map, indicating the presence of several in- 
dividual clouds or clumps with angular sizes of a few arc 
minutes which are unresolved by the telescope beam. We 
could confirm this clumpy structure in our follow-up syn- 
thesis observations of the area with the WSRT (Westmeier 
et al. 2005). Those clouds which could be spatially and kine- 
matically separated in our Effelsberg survey were included 
in the list of individual M31 HVCs in Table 2. The properties 
of the extra-planar gas and its possible origin in connection 
with the giant stellar stream (Ibata et al. 2001) of M31 are 
discussed in detail by Westmeier et al. (2005). 



6 COMPARISON WITH MODELS AND 
THEORETICAL PREDICTIONS 

One of the most fundamental questions arising from our ob- 
servations is how the detected HVC population around M31 

relates to different models derived from other observations 
or theoretical considerations. To investigate this question we 
will compare the observational parameters of the HVCs in 
a statistical approach with different models and theoretical 
predictions. 

6.1 The Local Group model of de Heij et al. 

(2002a) 

Based on the LDS and HIPASS, de Heij et al. (2002a) com- 
piled an all-sky catalogue of CHVCs. To study the distribu- 
tion and origin of the CHVC population, they constructed 
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Figure 11. Hi column density map of the extra-planax gas found near the eastern and south-eastern edge of the disc of M31. The 
centre of M31 is located at offset (0,0). Contours are drawn at 0.1, 0.25, 0.5, 0.75, 1, 1.25, 2.5, 5, 10, 20, 30, 40, 50, 60, 70, 80, and 
100 X 10^^ cm~^. Two example spectra at the positions labelled with A and B are shown on the right. The complex line profile at 
position A indicates the presence of compact substructures unresolved with the 9 arcmin HPBW of the Effelsberg telescope. The cloud 
at position B shows high intensities and narrow lines of only 15 kms'^ FWHM. 



several Local Group population models with different pa- 
rameters and compared the results with the observed pa- 
rameters of the CHVCs. These models were motivated by 
the suggestion of Blitz et al. (1999) that HVCs could be the 
gaseous counterparts of the 'missing' dark-matter satellites 
predicted by CDM structure formation simulations (Klypin 
et al. 1999; Moore et al. 1999). The free parameters of the 
models by de Heij et al. (2002a) include the Gaussian dis- 
persion, (7, of the radial distribution of CHVCs around their 
host galaxy in the range of cr = 100 kpc ... 2 Mpc, the slope, 
P, of the H I mass function in the range of /? = —2.0 ... —1.0, 
and the highest allowed H I mass, Mmax, of the clouds in the 
range of log (Mmax/M©) = 6.0 . . . 9.0. Using the test and 
the Kolmogorov-Smirnov (KS) test, de Heij et al. (2002a) 
determined the parameter space providing the best repre- 
sentation of the observational parameters of CHVCs. One of 
the best-fitting models is their model #9 with a = 200 kpc, 
P = -1.7, and M„,ax = iC M©. 



6.1.1 Construction of a model HVC population 

In order to compare our results with the parameters of 

model #9 of de Heij et al. (2002a), we can construct a 
model with similar statistical parameters and then correct 
for the incompleteness issues of our Effelsberg observations. 
According to model #9, the total number of HVCs around 
M31 should be iV = 750. This assumes equipartition of the 
entire model population between M31 and the Milky Way, 
which are both believed to have similar total mass (Evans & 
Wilkinson 2000). The spatial distribution of the clouds has 
to be spherically-symmetric with a Gaussian radial decline. 
Following the method introduced by Box & Muller (1958), 
we can assign each cloud in our model a radial distance from 
M31 of 

r = U/ —2(7^ In(xi) cos(27ra;2) (5) 



with a being the radial dispersion of the spatial distribution 
of our HVC population. Furthermore, xi and X2 must be two 
independent uniform random deviates in the range of < 
Xi < 1. The remaining two spatial coordinates, the angles 
ip and have to be determined such that the directions of 
the HVCs, as seen from the centre of M31, are distributed 
uniformly. 

The next step is to assign each of the model clouds a 
radial velocity Vrad. Again, we can use the method of Box & 
Muller (1958) to determine radial velocities with a Gaussian 
distribution around the systemic velocity of M31 of vmsi 
—300 kms^^. Hence, 

frad = VM31 + \/-2ct2 In(xi) 008(2^2) (6) 

with (Ji, being the line-of-sight velocity dispersion of the en- 
tire population of CHVCs which is expected to be governed 
by the details of the mass model assumed for M31. For the 
sake of simplicity we can use the same value as observed 
for the population of globular clusters around M31, namely 
a„ ^ 150 kms"^ (Perrett et al. 2002). 

In the final step we have to assign each model cloud 
an H I mass such that the postulated mass function is pre- 
served. In the model of de Heij et al. (2002a) a power-law 
mass function with a slope of /3 = —1.7 is assumed. We can 
reproduce this mass function via 

MHi = M^i„xi/(''+i\ (7) 

where Mmin is the lower end of the mass function, and x 
must again be a uniform random deviate in the range of 
< a; < 1. 



6.1.2 Consideration of the survey limitation 

With the ingredients listed above we are able to construct 
a model population of HVCs around M31 with parameters 
similar to model #9 of de Heij et al. (2002a) . Before we can 
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Figure 12. Projected radial distribution (left) and Hi mass spectrum (right) of HVCs around M31 according to our EflFelsberg obser- 
vations (solid line) and our model (dotted/dashed lines) based on model #9 of de Heij et al. (2002a). In both diagrams the dotted line 
corresponds to the original model, whereas the dashed line represents the remaining model population after correcting for the different 
incompleteness issues of our survey. 



compare this model population with our Effelsberg obser- 
vations we have to account for the different observational 
hmitations of our survey. The most obvious limitation is the 
survey area for which we can account by simply removing 
all clouds from the model which, when projected on the sky, 
arc located outside the survey area. Next, we have to remove 
all clouds which suffer from blending with Hi emission of 
the Galactic disc, namely all clouds with radial velocities of 
Vrad > -140 kms"\ 

A correct treatment of the flux detection limit of our 
survey is more complicated. To assess whether a cloud of 
certain mass can be detected with our survey we have to 
find a conversion between the Hi mass and the observed 
flux of a gas cloud. According to Schneider (1997), the total 
observed flux, F, of an unresolved gas cloud of H I mass Mh i 
and distance d is given by 

:4.24xl03^fA^-^ (8) 



m Jy km s ^ 

The Effelsberg spectra, however, arc given in terms of bright- 
ness temperature. Flux density, S, and brightness tempera- 
ture, Tb, can be converted via 



Tb 



2kfl 



(9) 



with k being the Boltzmann constant and Q the solid angle 
of the telescope beam. For a constant wavelength of A = 
21.1 cm and a circular telescope beam of HPBW 6, Eq. 9 
can be simplified to 



K 



S 



" ^ \ arcsec / 



(10) 



mjy beam 

Assuming the source being unresolved, we can combine Eq. 8 
and 10 to obtain a relation between the H I column density 
and the H I mass of a cloud, namely 



^=4.68xl0-^r^ 
cm-'' Mq \ arcsec kpc J 



(11) 



where we also made use of the conversion between the inte- 
grated brightness temperature, f Tb dv, and the H I column 
density, A^hi, of a cloud under the fundamental assumption 
that the gas is optically thin. Assuming a Gaussian line pro- 
file with dispersion a, we can derive the H I column density 



sensitivity of our Effelsberg survey from the brightness tem- 
perature sensitivity via 



ciTb I expf-I^i^'ld. 



= Cl 



4 In 0.5 



2a2 



TbAv « ciTbAv, 



(12) 



where Tb denotes the peak brightness temperature and wo 
the radial velocity of the spectral line. Av ~ 2.355 a is the 
FWHM of the spectral line. If A'hi is given in cm~^ and 
TbAv in Kkms~^, the corresponding conversion factor is 
Cl = 1.823 X lO'^*, again assuming that the gas is optically 
thin. 

As we have seen in Section 4.3, the detection rate of 
clouds at the 3.5 a level is about 50 per cent. Below 2cr, 
the detection rate drops to zero, whereas above 5cr we are 
virtually complete (see Fig. 3). Thus, it is reasonable to con- 
struct a soft cut-off level with a linear increase in detection 
probability -P(Tb) such that 



P(Tb) = { 







for Tb < 2 o- 



if- I for 2(7 < Tb < 5cT 
1 forTB>5(7 



(13) 



By converting the values for Tb into the corresponding Hi 
masses (using Eq. 11 and 12) we are able to exclude clouds 
from the model with a certain probability expressed by 
Eq. 13 to account for the detection limit of our Effelsberg 
data. 



6.1.3 Comparison between model and data 

To minimize statistical errors we calculated the mean pa- 
rameters of 1000 model HVC populations with the speci- 
fications discussed above. Each model originally contained 
750 HVCs of which on average 595 were excluded from the 
model population by the incompleteness issues of our Effels- 
berg survey, leaviug us with about 155 expected HVC de- 
tections within our survey area. This number immediately 
illustrates that model #9 of de Heij et al. (2002a) overesti- 
mates the expected number of HVCs detected in our survey 
by about a factor 10. 
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Figure 13. Normalised cumulative distribution of HVCs as a function of projected distance, dproj • from the centre of M31 (left) and H I 
mass, Mhi (right). The solid line shows the 17 detected HVCs from our survey, whereas the dashed curve is the expected distribution 
according to our model based on model #9 of de Heij et al. (2002a). 



The resulting distribution of HVCs as a function of 
projected distance, dproj, from the centre of M31 is plot- 
ted in the left diagram of Fig. 12. The dotted curve repre- 
sents the original model containing 750 HVCs. The dashed 
curve shows the final distribution after accounting for the 
incompleteness effects. The distribution drops abruptly be- 
yond dproj ^ 150 kpc as a result of the limited survey re- 
gion. The observed radial distribution is plotted for com- 
parison as the solid Une. It is obvious that both the total 
number and the radial distribution of the observed HVCs 
around M31 deviate significantly from the model. The ob- 
served HVCs are much closer to M31 than expected from the 
model, whereas no HVCs are observed at larger projected 
distances where the model still predicts several detections. 
To quantify these differences in radial distribution we carried 
out the Kolmogorov-Smiriiov (KS) test on the normalised 
cumulative radial distribution of the modelled and observed 
HVCs as plotted in the left panel of Fig. 13. The resulting 
KS statistic for the two distributions \s D = 0.32, hence 
rejecting the hypothesis that both samples are drawn from 
the same radial distribution on the 95 per cent confidence 
level. 

The deviations are similar when considering the Hi 
mass spectrum. The modelled mass spectrum is shown in 
the right panel of Fig. 12. Again, the dotted curve is the 
original mass spectrum, whereas the dashed curve has been 
corrected for observational constraints. The modelled distri- 
bution essentially follows the power law defined by de Heij 
et al. (2002a). In contrast, the observed Hi mass spectrum 
(solid line) does not appear to follow a power law at all. In- 
stead, it peaks at Mhi « 5 x 10^ Mq, indicating a lack of 
low-mass clouds compared to the model. To compare the ob- 
served and modelled mass spectra we again used the KS test 
on the normalised cumulative mass distributions as plotted 
in the right panel of Fig. 13. Again, the resulting KS statis- 
tic is Z) = 0.32, thus rejecting the modelled mass function 
on the 95 per cent confidence level. 

Obviously, the model proposed by de Heij et al. (2002a) 
does not provide a good and accurate description of our 
data. The model was motivated phcnomenologically based 
on the population of CHVCs observed all over the sky. How- 
ever, both the total number and the radial distribution of the 
HVCs detected around M31 do not agree with the predic- 



tions of the model. This result confirms our conclusion that 
HVCs/CHVCs arc much closer to their host galaxies than 
previously assumed, implying not only smaller projected dis- 
tances but also lower H I masses. Hence, a large fraction of 
the HVC population around M31 may have Hi masses be- 
low our 3 a detection limit of 8 x 10* Mq and could have 
remained undetected. 



6.1.4. Modification of the original model 

To study the radial distribution of the HVCs around M31 
in more detail, we modified the model of de Heij et al. 
(2002a) in order to determine the radial scale length, a, of 
the modelled HVC population which provides the best fit 
to the observed distribution of projected radii. In practice, 
we calculated several models covering the entire range of 
(T = 20 . . . 200 kpc, but leaving all other parameters (total 
number of clouds. Hi mass spectrum, etc.) fixed. Again, we 
calculated 1000 individual model populations for each value 
of cr to keep statistical errors small. We then compared the 
projected radial distribution of the model clouds with the 
observed one, again using the KS test. The results are pre- 
sented in Fig. 14. 

The left diagram of Fig. 14 shows the KS probability for 
the hypothesis that the projected radial distributions of ob- 
served and modelled HVCs agree as a function of radial scale 
length, (7, of the model population. The probability curve 
shows a very pronounced and narrow peak at cr w 49 kpc 
where the probability exceeds 80 per cent. Towards larger, 
and in particular towards smaller, scale lengths the KS prob- 
ability drops very quickly. For o = 100 kpc the probability 
is of the order of 15 per cent, and it decreases to only about 
5 per cent for the original de Heij et al. (2002a) model of 
(T = 200 kpc. The normalised cumulative radial distributions 
of model and data for the best-fitting case, a = 49 kpc, are 
plotted in the right panel of Fig. 14. In comparison to the 
original model of de Heij et al. (2002a) in Fig. 13 the ra- 
dial distribution is described much better, although model 
and data do not match in every detail. The observed HVC 
population shows a significantly steeper rise in the radial 
profile than the model population. This steep rise is caused 
by the fact that 6 HVCs belong to the complex of clouds 
near the south-eastern edge of M31, resulting in similar pro- 
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Figure 14. Left: Kolmogorov— Smirnov probability for the hypothesis that the observed radial distribution of HVCs around M31 follows 

a Gaussian distribution with radial scale length a. There is a clear peak at cr = 49 kpc with a probability in excess of 80 per cent. 
Right: Comparison of the normalised cumulative distributions of observed (solid line) and modelled (dashed line) HVCs as a function of 
projected distance from the centre of M31 for the best-fitting radial scale length of cr = 49 kpc. The KS test in this case yields a statistic 
of -D = 0.15, resulting in a probability of p = 81 per cent that both samples were drawn from the same distribution. 



jected radii. The appearance of this complex suggests that 
its constituents have a common origin and should not be 
treated as separate objects. 

In summary, the original Local Group population model 
of de Heij et al. (2002a) can be excluded with high confi- 
dence. Instead, our results support a more compact popu- 
lation with a Gaussian radial scale length of cr ~ 50 kpc. 
However, the radial scale length is not strictly confined by 
our observations. Even for a = 100 kpc we still obtain KS 
probabilities of about 15 per cent for agreement of our ob- 
servations with the model. 

6.2 Comparison with the satellite galaxies of M31 

Like the Milky Way, M31 hosts an extended group of satel- 
lite galaxies of which M33 is the most massive one. An 
overview of our current census of the M31 group is given by 
McConnachie & Irwin (2006) and Koch & Grebel (2006). 
In addition, several very faint dwarf galaxies, named An- 
dromeda IX-XIII, have recently been discovered by Zucker 
et al. (2004, 2006) and Martin et al. (2006). These detections 
have shown that a large number of low-mass satellites with 
very low surface brightness may still be awaiting discovery 
and that the M31 group, as presently known, is presumably 
incomplete. 

The radial distribution of satellite galaxies around M31 
has been reasonably well determined (see, e.g., Koch & 
Grebel 2006 for an overview) , although the relative errors of 
galactocentric distances can be significant for the innermost 
satellites of M31. Due to the lack of distance information, 
the spatial distribution of the HVCs around M31, however, 
is completely unknown. Hence, the only feasible approach 
is to compare the projected positions of HVCs and satellite 
galax;ies by discarding the distance information available for 
the M31 satellites. A similar problem arises when trying 
to compare the mass spectrum of both populations, as for 
most of the M31 satellite galaxies dynamical studies are not 
available. There have been Hi studies, however, for several 
gala^xies, although in four cases only upper limits resulting 
from non-detections are available (see Westmeier 2007 for a 
complete list of Hi observations and references). 



In Fig. 15 the logarithm of the Hi mass, Mni, of HVCs 
(open squares) and satellite galaxies (filled squares) is plot- 
ted versus the logarithm of the projected distance, dproj • 
Only those satellite galaxies within a projected radius of 
300 kpc around M31 have been considered as they axe ap- 
proximately located within the expected virial radius of M31 
(see, e.g., Kravtsov et al. 2004). There appears to be a dis- 
tinction between HVCs and satellite gala^xies. The HVCs 
are concentrated in the lower left area of the diagram with 
lower Hi masses and smaller projected distances, whereas 
the galaxies tend to concentrate in the upper right region 
of the diagram with typically higher Hi masses and larger 
distances. Apparently, HVCs and satellite galaxies have a 
completely different radial distribution. All HVCs are found 
within a projected radius of about 50 kpc, whereas the satel- 
lite galaxies of M31 are spread across the entire range of dis- 
tances out to 300 kpc. The Kolmogorov-Smirnov test yields 
a statistic oi D = 0.78, resulting in a very low probability of 
only 1.6 X 10~^ that HVCs and satellite galaxies are drawn 
from the same projected radial distribution. 

When interpreting the mass-distance diagram we have 
to keep in mind that the incompleteness and selection effects 
discussed above apply. As a consequence, the separation be- 
tween HVCs and galajcies in Fig. 15 may be an artefact of 
incomplete data. For several satellites of the M31 group, H I 
measurements are either not available or only upper mass 
limits have been derived from non-detections. As this will 
particularly affect those galaxies with intrinsically low Hi 
masses, a large number of data points could be missing in 
the low-mass domain, and there could be a more signifi- 
cant overlap between galaxies and HVCs than suggested by 
Fig. 15. 

If we assume the difference between HVCs and galax- 
ies in Fig. 15 to be real, what could the physical reason for 
the observed discrepancies be? If both HVCs and satellite 
galaxies were the visible counterparts of dark matter satel- 
lite haloes, the differences in their radial distribution would 
be surprising. In this ease, a physical mechanism would be 
required to prevent the presence of dark-matter-dominated 
neutral gas clouds at large galactocentric distances from 
M31, e.g. through ionisation of the gas. This scenario will 
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Figure 15. Hi mass versus projected distance of HVCs (open 
squares) and satellite galaxies (filled squares) around M31. Al- 
though both groups tend to populate different regions of the dia- 
gram, it is important to note that for only a fraction of the M31 
satellites Hi studies are available. 

be discussed in more detail in the following Sect. 6.3. An- 
other reason for the discrepancies in the radial distribution 
of HVCs and satellite galaxies could be a different origin of 
both groups of objects. If HVCs were primarily of tidal ori- 
gin or the result of condensing halo gas, gaJactic winds, or 
outflows, their proximity to M31 would be a natural conse- 
quence. 

6.3 Comparison with CDM simulations 

Kravtsov et al. (2004) carried out structure formation simu- 
lations of Galaxy-sized dark-matter haloes in the framework 
of ACDM cosmology. They also included a simple model 
of star formation and gas dynamics in their simulations. A 
comparison with their results reveals significant differences 
in the radial distribution of both HVCs and satellite galax- 
ies compared to the predicted population of dark matter 
satellites around Galaxy-sized haloes. Within a radius of 
200 kpc the mean distance of dark matter satellites of 
about 120 kpc is noticeably larger than that of the M31 
satellite galaxies of 85 kpc and significantly larger (as- 
suming h 0.7) than the maximum projected distances of 
about 50 kpc observed for HVCs. Kravtsov et al. (2004) can 
partly solve this problem by including a simple model of 
star formation in dark matter haloes in their simulations. 
According to their results only a small fraction of satellites 
have actually formed stars to become luminous galaxies, 
whereas most haloes remained dark. Their model can not 
only reproduce the radial distribution and number of satel- 
lite galaxies observed around the Milky Way, but also the 
circular velocity function and the observed morphological 
segregation. 

To investigate whether the remaining dark satellites in 
their simulations can be identified with the population of 
HVCs around M31 or the Milky Way, Kravtsov et al. (2004) 
also determined the gas mass associated with each dark mat- 
ter halo in their simulation. Within a radius of 200 kpc 
they find a total gas mass of Mg « 2 x 10^ M0, resulting 
in a corresponding neutral gas mass of Mhi w 2 x 10* M© 
if a neutral gas fraction of 10 per cent is assumed (e.g.. 



Maloncy & Putman 2003). This value is one order of mag- 
nitude higher than the total H I mass of about 2 x 10^ Mq 
of the 17 HVCs identified in our Effelsberg survey, including 
Davies' Cloud. We should note, however, that some more 
neutral gas is present in the form of diffuse, extra-planar 
gas which is not resolved into individual clouds with our 
9 arcmin HPBW. 

The total number of gaseous dark haloes with Mg > 
lO** M0 (corresponding to Mhi > 10^ Mq if a neutral hy- 
drogen fraction of / = 0.1 is assumed) in the simulations of 
Kravtsov et al. (2004) is about 50 . . . 100, but for the cen- 
tral 50 kpc their models predict only 2 ... 5 such clouds. 
They conclude that part of the HVCs found around M31 by 
Thilker et al. (2004) could be primordial dark matter haloes, 
whereas others could instead be of tidal origin similar to the 
Magellanic Stream around the Milky Way. In our Effelsberg 
survey, however, we have not detected any additional HVCs 
beyond a projected radius of about 50 kpc with a 3a de- 
tection limit of 8 X lO'* M0. This discrepancy is surprising, 
given that more than 90 per cent of the predicted gaseous 
haloes should be located outside a radius of 50 kpc. 

Because of the limited azimuthal range of our survey 
region we cannot assume to have detected all of the ex- 
pected gaseous haloes. At a projected distance of 50 kpc 
our local completeness function drops below 50 per cent (see 
Fig. 2), but none the less we still cover more than one third of 
the total area within a projected radius of 100 kpc. Within 
a radius of 100 kpc Kravtsov et al. (2004) expect about 
35 per cent of all dark matter haloes, resulting in 18 ... 35 
expected gaseous haloes with Mg > 10® Mq , about a third of 
which should be located within our survey area, resulting in 
a total number of expected HVCs in our survey in the range 
of about 6 ... 12 (not taking into account the details of the 
radial distribution of dark matter haloes). This is less than 
the 17 HVCs identified in our data, supporting the concept 
of Kravtsov et al. (2004) that some of the detected HVCs 
may represent tidally stripped gas originating from satellite 
galaxies of M31. However, we would expect the gaseous pri- 
mordial dark matter haloes to be spread across the entire 
survey area, whereas the detected HVCs are highly concen- 
trated in the central region around M31. 

One possible solution to this problem could be ionisa- 
tion. With their masses and sizes, the HVCs observed around 
M31 are similar to the HVCs considered by Sternberg et al. 
(2002) in their hydrostatic simulations of compact HVCs. 
For their circumgalactic model, Sternberg et al. (2002) as- 
sumed spherical, pressure-confined, dark-matter-dominated 
clouds at typical distances from the Galaxy of the order of 
150 kpc. The original aim of the simulations was to explain 
the population of CHVCs as defined by Braun & Burton 
(1999) and to understand their observational parameters as 
determined by Braun & Burton (2000) and Burton et al. 
(2001). For the circumgalactic case, Sternberg et al. (2002) 
predict typical Hi masses of 3 x 10^ M© and a radial scale 
length of the H I column density distribution of the clouds of 
0.5 kpc. The typical peak Hi column density of their model 
clouds is 5 X lO'^^ cm~^. These parameters are in very good 
agreement with those of the HVCs around M31 as observed 
with the Effelsberg telescope and the WSRT (Westmeier 
et al. 2005). 

The HVCs modelled by Sternberg et al. (2002) are grav- 
itationally dominated by a dark matter halo for which a 
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Burkert profile (Burkert 1995) was assumed. For the circum- 
galactic case, Sternberg et al. (2002) assumed a virial mass 
of the halo of the order of 10* Mq. This implies a virial-to- 
Hi meiss ratio of the order of a = 300, which corresponds 
to a very small neutral gas fraction of only 0.3 per cent. 
Despite the substantial dark matter content, the halo would 
not be massive enough to retain the warm gas. Therefore, 
the clouds have to be additionally stabilised by the pressure 
of an external medium of at least P/k « 50 K cm~^. The ex- 
pected circumgalactic hot ionised corona around the Milky 
Way and M31 could provide the required external pressure 
stabilisation. In this case, the central gas density of the mod- 
elled HVCs would be A^h « 2 x 10^^ cm^'^ which is again 
in excellent agreement with the observed mean neutral gas 
densities of the HVCs around M31 (Westmeier et al. 2005). 

This remarkable agroomont between the hydrostatic 
simulations of dark-matter-dominated HVCs by Sternberg 
et al. (2002) and the observed parameters of the HVCs de- 
tected in our WSRT observations suggests that some of the 
clouds could indeed be primordial, dark-matter-dominated 
clouds as originally assumed by Blitz et al. (1999) and Braun 
& Burton (1999). This could particularly be true for those 
clouds at larger projected distances from M31 which are 
completely isolated from all other HVCs and from M31 and 
any of its known satellite galaxies. The required pressure 
confinement by the corona of M31 could also explain the 
head-tail morphologies observed for some of these clouds by 
Westmeier et al. (2005). As the clouds are moving with high 
velocities through their environment their outer parts will 
be distorted and stripped by the ram-pressure of the ambi- 
ent medium (for hydrodynamic simulations see, e.g., Quilis 
& Moore 2001; Konz et al. 2002). 

The simulations of Sternberg et al. (2002) now allow us 
to find a plausible explanation for the concentration of HVCs 
close to M31 in contrast to the simulations of Kravtsov 
et al. (2004). At larger distances from M31 the density of 
the circumgalactic environment is expected to decrease. At 
some distance, the density could drop below the value of 
P/k « 50 K cm~"^ which is required in addition to the gravi- 
tational potential to retain the warm gas of the HVCs. This 
would have two consequences: First, the clouds would not be 
stable any more because the gas would no longer be bound. 
Second, the gas density would become too low to provide suf- 
ficient shielding against the ionising intergalactic radiation 
field. Consequently, the gas would be mainly ionised and no 
longer detectable in the 21-cm emission of neutral hydrogen. 
If this scenario is correct, hundreds of mainly ionised or pure 
dark-matter satellites should be present around large spiral 
galaxies, such as the Milky Way and M31, which would not 
be detectable in the 21-cm line of neutral hydrogen. 

The dark matter mini-halo model for HVCs by Stern- 
berg et al. (2002) provides a plausible explanation for some 
of the HVCs near M31, in particular those at larger pro- 
jected distances which are isolated in position- velocity space 
from both the stellar stream and any known dwarf compan- 
ion of M31. Tidal interaction during their evolution as well 
as instability and ionisation at larger distances from M31 
can explain the small number of HVCs and their proximity 
to M31. At the same time, HVCs would be sensitive test par- 
ticles to probe the physical conditions in the circumgalactic 
environment of large spiral galaxies. 



7 ORIGIN OF THE HVCS NEAR M31 

As discussed before, our observations seem to favour two 
possible scenarios for the origin of the HVCs found near 
M31: They could be tidally stripped gas clouds or primor- 
dial, dark-matter-dominated satellites. In this section we will 
discuss the evidence for these two scenarios and a few other 
popular hypotheses on the origin of HVCs. 

Several of the clouds could be the result of tidal or ram- 
pressure stripping in connection with the interaction be- 
tween M31 and its present or former satellite galaxies. This 
could be the case for the concentration of HVCs (number 1- 
6 in our catalogue) detected near the south-eastern edge of 
the H I disc of M31 which spatially and kinematically over- 
lap with the giant stellar stream (Ibata et al. 2001; Ferguson 
et al. 2006). The nature of this HVC complex was discussed 
in detail by Westmeier et al. (2005) based on high-resolution 
WSRT observations. The HVCs are characterised by a par- 
ticularly complex velocity structure (see Fig. 11) suggestive 
of turbulent motions. Unfortunately, recent numerical sim- 
ulations of the stellar stream (Fardal et al. 2006, 2007; Mori 
& Rich 2008) only focussed on non-coUisional particles and 
did not include gas, so that an association of the HVCs with 
the stellar stream can not be confirmed independently. An- 
other specific candidate for tidal or ram-pressure stripping is 
HVC number 14 which is located close to NGC 205. Again, 
the radial velocity of the HVC is similar to the velocities of 
the H I gas detected in NGC 205 (see Westmeier et al. 2005 
for details). 

A tidal origin would naturally explain the observed con- 
centration of HVCs in the vicinity of M31 where tidal forces 
and halo gas density are increased. Another consequence of 
the tidal scenario would be moderately high motallicities 
due to processing of the HVC gas in satellite galEixies and 
possible mixing with enriched gas from the halo or disc of 
M31. The only confirmed example of tidally stripped high- 
velocity gas in the vicinity of the Milky Way is the Magel- 
lanic Stream (including the Loading Arm) which shows typi- 
cal metallicities of about 0.1 .. . 0.4 (e.g., Wakker 2001; Sem- 
bach et al. 2001). Other large HVC complexes near the Milky 
Way show comparable metallicities of the order of 0.1 (see 
Wakker 2001 for a comprehensive review), demonstrating 
that the gas must have been processed and moderately en- 
riched in the past. Direct metallicity measurements of HVCs 
near M31, however, will be hampered by the lack of suitable 
background sources for absorption line spectroscopy. 

As discussed in Sect. 6.3, some of the HVCs around 
M31 could also be primordial dark- matter haloes left over 
from structure formation in the Local Group. This could 
in particular be the case for isolated HVCs at larger pro- 
jected distances from M31. In this case, we would expect 
to observe primordial metal abundances unless the gas has 
been enriched through interaction with halo gas or luminous 
dark-matter satellites. As mentioned before, another conse- 
quence would be the presence of ionised or pure dark-matter 
satellites out to distances of several hundred kpc from M31 
or the Milky Way. Depending on the nature of dark matter, 
pure dark-matter haloes would be very hard to detect, but 
they could leave noticeable imprints (e.g., Wang et al. 2007). 

The discovery of HVCs and several new satellite galax- 
ies near M31, however, cannot directly solve two fundamen- 
tal problems in connection with CDM structure formation 
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scenarios. First of all, the number of gaseous and luminous 
satellites would still be one order of magnitude below the 
predictions of CDM simulations, and we would still have to 
assume the additional presence of hundreds of pure dark- 
matter haloes to solve the 'missing satellites' problem. Sec- 
ond, The discovery of new satellites does not solve the prob- 
lem that satellite galaxies around the Milky Way and M31 
appear to be arranged in great planes with large inclination 
angles (e.g., Kroupa et al. 2005; Koch & Grebel 2006). Even 
if the new satellites turned out to have a spherical distribu- 
tion, the anisotropy of the most massive satellites would re- 
main. This anisotropy has been explained with the preferred 
infall of dark-matter haloes along the filaments of the cosmic 
web and subsequent evolution of the satellite orbits within 
the triaxial potential of their host galaxy's dark-matter halo 
(Libeskind et al. 2005; Zentner et al. 2005). An alternative 
hypothesis claims that many satellite galaxies around the 
Milky Way and M31 arc tidal dwarf galaxies instead of cos- 
mological dark- matter satellites (Metz et al. 2007). 

Another possible explanation for the more isolated 
HVCs found around M31 is a tidal origin in connection with 
a previous close encounter between M31 and M33. Accord- 
ing to recent test particle simulations by Bekki (2008) , tidal 
interaction between the two galaxies during their first en- 
counter between 4 and 8 Ga ago would have stripped gas 
from M33 which could account for the Hi bridge found 
by Braun & Thilker (2004) between M31 and M33. Bekki 
(2008) suggests that tidally stripped gas from M33 could also 
have formed some of the HVCs near M31. Consequently, the 
HVCs would have metallicities similar to those observed in 
the outer gas disk of M33. The limitations of current obser- 
vations and simulations, however, do not allow us to further 
investigate this interesting option. 

Another hypothesis suggests that HVCs are the result of 
galactic outflows due to individual supernovae or the galac- 
tic fountain process (e.g., Bregman 1980; Booth & Theuns 
2007). This hypothesis appears to be inconsistent with the 
typically low metallicities of Galactic HVCs of the order of 
10 per cent of the solar value (see Wakker 2001 for a com- 
prehensive review). However, the measured values vary over 
some range, indicating that the HVCs around the Milky 
Way do not form a homogeneous population. In addition, 
metallicities have only been measured for the large HVC 
complexes near the Milky Way, whereas the metallicities 
of most CHVCs are unknown. At the same time, we failed 
to detect the expected large population of CHVCs around 
M31 in our Effelsberg survey, suggesting that CHVCs could 
be intrinsically compact objects at small distances from the 
Milky Way and M31. Hence, CHVCs could be promising 
candidates for gas ejected by galactic outflows. 



8 SUMMARY AND CONCLUSIONS 

We used the 100-m radio telescope at Effelsberg to map a 
large area around the Andromeda Galaxy, M31, in the 21- 
cm line emission of neutral atomic hydrogen. Our survey 
extends out to a projected distance of about 140 kpc in the 
south-eastern direction (equivalent to about two thirds of 
the projected distance towards M33) and about 70 kpc in 
the north-western direction. With this map layout we are 
able to flu the previously existing gap between the outer 



boundary of the GBT survey of Thilker et al. (2004) at about 
50 kpc projected distance from M31 and the upper limit of 
about 150 kpc for the distance of HVCs as derived from 
the non-detections in nearby galaxy groups (Zwaan 2001; 
Braun & Burton 2001; Pisano et al. 2004). The achieved 
spectral baseline RMS is 45 mK at 2.6 kms~^ velocity res- 
olution, corresponding to a 3 cr H I column density detec- 
tion limit of 2.2 x 10^* cm~^ for the warm neutral medium 
(Av = 25 kms"'^ FWHM). This translates into an Hi mass 
sensitivity of 8 x 10'' M©. 

In total, we detected 17 individual HVCs and several 
regions of more extended extra-planar gas all around M31. 
The discrete clouds are predominantly unresolved by the 
HPBW of 9 arcmin and characterised by typical H I masses 
of a few times 10 ' Mg . We did not detect any clouds be- 
yond a projected distance of about 50 kpc, suggesting that 
HVCs are generally found in proximity of their host galax- 
ies. In particular, we did not find an extended populations 
of hundreds of CHVCs as observed around the Milky Way, 
suggesting that the Galactic CHVCs are intrinsically small 
clouds in the immediate vicinity of the Milky Way. 

A comparison with the Local Group population model 
of CHVCs, as proposed by de Heij et al. (2002a), reveals 
that their best-fitting model #9 is discarded by our data 
with high confldence. Neither the observed projected radial 
distribution nor the H I mass function of HVCs around M31 
can be explained by the model. Instead, we find that a Gaus- 
sian radial scale length of the order of 50 kpc can best ex- 
plain the observed projected distribution of HVCs around 
M31. In addition, we find that the HVCs are also distinct 
from the M31 satellite galaxies through typically lower Hi 
masses and smaller projected distances from M31. 

CDM-bascd structure formation simulations by 
Kravtsov et al. (2004) suggest that about 50 to 100 dark- 
matter haloes with total gas masses of greater than 10® M© 
should exist within about 300 kpc of M31. Only 2 to 5 
of these haloes should be located within 50 kpc of M31, 
suggesting that some of the HVCs found near M31 could 
instead be tidally stripped gas from present or former 
satellite galaxies of M31. This idea is supported by our 
high-resolution follow-up observations of several HVCs 
(Westmeier et al. 2005). 

The lack of detections beyond 50 kpc projected rar 
dius, however, is in conflict with the predictions made by 
Kravtsov et al. (2004). A possible explanation of this dis- 
crepancy could be ionisation as a result of decreasing pres- 
sure of the ambient coronal medium at larger distances from 
M31, as suggested by hydrostatic simulations of Sternberg 
et al. (2002). An important consequence of this scenario 
would be the presence of hundreds of mainly ionised or pure 
dark-matter satellites near large spiral galaxies, such as the 
Milky Way and M31, which would be undetectable in the 
21-cm line of neutral hydrogen. Finding these almost invisi- 
ble satellites would be an important observational result in 
support of CDM cosmologies. 

Another promising scenario recently discussed by Bekki 
(2008) is a possible tidal interaction between M31 and M33 
during their previous encounter. This would have stripped 
gas from M33 which could have formed some of the HVCs 
observed near M31. More detailed and extensive observa- 
tions and simulations will be required to further investigate 
this interesting scenario. 
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With our Effelsberg Hi survey of M31 we have shown 
that HVCs are most hkely concentrated around the large 
spiral galaxies with typical distances of no more than a few 
10 kpc. It is likely that different physical processes, such as 
tidal stripping, accretion of primordial dark-matter haloes, 
or galactic outflows, have contributed to the HVC popula- 
tions of M31 and the Milky Way. The general volatility of 
these processes could naturally explain some of the differ- 
ences between the HVC populations of M31 and the Milky 
Way, for example the extended and complex gaseous streams 
of the Magellanic Clouds which have no counterpart in M31. 
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